Closed-loop Quantum Parameter Estimation: Spins in a Magnetic Field 
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We present an experimental demonstration of closed-loop quantum parameter estimation in which 
real-time feedback is used to achieve robustness to modeling uncertainty. By performing broadband 
estimation of a magnetic field acting on hyperfine spins in a cold atom ensemble, we show that 
accuracy is not compromised by fluctuations in total atom number even though the measured 
signal in our canonical configuration depends only on the product of the field and atom number. 
This methodology could be essential for efforts to utilize conditional squeezing in spin-resonance 
measurements. 
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PACS numbers: 06.20.-f, 32.80.Pj, 32.80.Qk 

Optimal design of experimental procedures and data 
analysis strategies can often be accomplished using tech- 
niques from parameter estimation theory. Such an ap- 
proach can be essential for obtaining acceptable perfor- 
mance in scenarios where modeling is subject to some 
degree of uncertainty, and has seen widespread use in 
fields that lie near the interface of physics and informa- 
tion science. Prominent current examples of such fields 
include metrology, optical communication, and computa- 
tion with novel substrates. As micro- and nano-scale sys- 
tems with manifestly non-classical behavior have gained 
importance in these areas, researchers have devoted in- 
creasing attention to the extension of parameter estima- 
tion methodologies to problems involving quantum me- 
chanical inference and dynamics. Our aim in this ar- 
ticle is to establish that real-time feedback plays a cen- 
tral role in the quantum regime — as it does in classical 
scenarios — enabling robust parameter estimation in the 
presence of significant modeling uncertainty. 

In quantum parameter estimation the central objec- 
tive is to extract information about a static or time- 
dependent parameter in a Hamiltonian, via direct or in- 
direct measurements performed on a system that evolves 
according to this Hamiltonian 0, 0, • An elementary 
example of such a process is estimating the amplitude of a 
magnetic field, b(t), by observing short-time Larmor pre- 
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FIG. 1: Quantum parameter estimation of a magnetic field 
by observing ensemble spin dynamics according to both an 
open (A) and closed (B) loop experimental methodology. 



cession of a spin ensemble ^ (here we will consider an 
ensemble consisting of hyperfine spins in a cloud of laser- 
cooled atoms). The field-spin interaction is described by 
a magnetic dipole Hamiltonian, 



H[t;b(t)] = hjh(t) ■ F 



(1) 



where 7 is the gyromagnetic ratio and F is the total an- 
gular momentum operator for the ensemble. 

A canonical estimation procedure is depicted in Fig. 
mA). An atomic spin ensemble is prepared such that 
its net magnetization, or Bloch vector, (F), achieves a 
coherent spin state along the a;-axis. In the presence 
of an external magnetic field, b(i), the bulk magnetiza- 
tion precesses around b with instantaneous frequency, 
= 7|b(t)|. Generally, it is arranged such that 
b(t) = b(t)y lies along the y-axis so that the spin state 
can be resolved from a continuous measurement of (F z ) 
5, 6]. An estimate, b, of the external magnetic field is 
then obtained (by regression or filtering) from the small- 
angle relation 
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(2) 



This procedure can achieve Heisenberg limited sensitivity 
by exploiting conditional spin-squeezing 0, @ . The 
problem with this approach is that it requires accurate 
knowledge of the net magnetization, |F(t)|, a quantity 
that unfortunately varies due to decoherence and shot- 
to-shot fluctuations in the atom number. Uncertainty in 
|F(i)| directly translates into uncertainty in b(t). 

Fig. ^B) depicts an alternative estimation procedure 
that is robust to fluctuations in |F(t)| 0. The spin en- 
semble and (F z ) measurement are situated within a feed- 
back control loop that attempts to stabilize (F z ) to a ref- 
erence value, r(t) = 0. In the presence of a time-varying 
magnetic field signal, h(t) — b(t)y, the controller imposes 
a compensating field, b c (t) ~ — b(t)y, to try to suppress 
the atomic Larmor precession. The closed-loop estimate 
is then given by b(t) = — b c , and its accuracy is deter- 
mined by the controller's ability to respond promptly and 
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FIG. 2: Schematic of our apparatus for closed-loop parameter 
estimation using a cold atom ensemble to determine h(t) in a 
magnetic dipole Hamiltonian, H(t) = ?ryb(t) • F. 



accurately to changes in b(i). In this work we demon- 
strate that standard design techniques enable the imple- 
mentation of feedback controllers with excellent tracking 
and high robustness to fluctuations in |F(i)|. This il- 
lustrates new utility for real-time feedback in cold atom 
physics, as previous investigations have focused on active 
control of atomic motion [llj, |l2j • 

Fig. |21 provides a schematic of our experimental ap- 
paratus implemented according to the design in Fig. 
H^B). It consists of a cold atom ensemble, a Faraday 
polarimeter for continuously probing the atomic Bloch 
vector, and a high-speed Helmholtz coil along the y-axis 
for applying feedback. The spin system is provided by 
the 6 2 S!/2(F=4) ground state hyperfine manifold in Cs, 
which contains (2F + 1) = 9 energetically degenerate 
Zeeman sublevels in zero field. Therefore, the net an- 
gular momentum of the polarized spin state is given by, 
Nhy/F(F + 1) w ANU, for an ensemble of N atoms. 

Cold atom samples were obtained by loading N ~ 10 9 
neutral 133 Cs atoms into a magneto-optical trap (MOT) 
from a ~ 1 x 10~ 8 Torr background vapor, using optical 
trapping beams (30 mW each, 2.5 cm diameter) derived 
from an injection-locked diode laser. The atoms were 
cooled to 1 fiK via a 5 ms cr+ /<r_ polarization gradient 
cooling phase, and a coherent spin state was produced 
by optical pumping with a + polarized light (100 fxW 
with a 6.2 mm Gaussian waist) tuned to the 6 2 S!/ 2 (F=4) 
— >6 2 P 3 /2(F'=4) hyperfine transition. A 45 mW, 2.5 cm 
diameter re-pumping beam was used. 

Continuous weak measurement of F z was implemented 
using a free-running diode laser, blue-detuned from the 
(F=4)->(F'=5) transition by A=2 GHz. The 65 
probe beam was linearly polarized by a high extinc- 
tion (>10 6 ) Glan-Thompson prism polarizer and mode- 
matched to the Gaussian waist of the atomic spatial dis- 



tribution. Relaxation of the coherent spin state due to 
the probe light was measured to be T% = 11.2 ms. 

Faraday rotation of the probe light was detected using 
a balanced polarimeter (1 MHz detector bandwidth) con- 
structed from a high extinction Glan-Thompson polariz- 
ing beam splitter. This configuration produces a pho- 
tocurrent, y(t), proportional to the ^-component of the 
spin angular momentum 0, 0, llOj , 



(3) 



y(t)=2VM{F z (t))+C(t) 



where M is the measurement strength and reflects 
measurement noise. Background magnetic fields were 
nulled to <100 /iG with large (d = 1 m) external three- 
axis Helmholtz coils by Larmor precession measurements. 

The feedback system in Fig. [21 utilizes the photocur- 
rent, y(t), to control the strength of a Hamiltonian, 



H c (t) = h-yb c [y(t)]F y 



(4) 



that rotates the Bloch vector around the y-axis. It is the 
job of the controller to determine the appropriate feed- 
back strength, b c (t) — flu(t), based on the observation 
y(t). For a linear controller, 



u{t)= / C(t- T )y(T)dT, 
Jo 



(5) 



where C(t) must be designed to satisfy tracking and 
robustness objectives. The controller output, u(i), is 
used to program a current supply that drives a y-axis 
Helmholtz coil surrounding the atomic sample. This re- 
quires accurate calibration (obtained in our case from 
Larmor frequency measurements) of the gain, 0, from 
current supply programming voltage, u(t), to the feed- 
back field, b c (t) . 

In control theory, it is customary to design C (t) in the 
frequency rather than time domain by taking Laplace 
transforms of the functions, C(t) — > C(s), y(t) — > y(s), 
etc., where s = ju> and j = This results in the 

closed-loop frequency response of the control system, 



T(s) 



y( s ) 



C(s)P(s) 



r(s) l + C{s)P(s) 



(6) 



where P(s) = y(s)/u(s) is the transfer function from the 
feedback coil programming voltage to the photocurrent. 
The tracking objective is to adjust the gain and phase of 
C(s) such that T(s) ~ 1 over as large a frequency range 
as possible. Standard results from robust loop-shaping 
theory 0,0] lead to a controller design, 



coo 



Q(s) 



l + P(s)Q(s)' 



Q = 



W(s) 

Pm P (s) 



(7) 



where W(s) is a stable, strictly proper [W(s — > oo) = 0] 
function that adjusts the bandwidth of C(s) and P mp (s) 
is the minimum phase contribution to P(s). 
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FIG. 3: (A) Measured (points) and fitted (solid line) plant transfer function, P(s), between the feedback signal, u(t), and 
polarimeter photocurrent, y(t), as well as the designed controller, C(s). (B) Measured (points) and calculated (solid line) 
closed-loop transfer function, T(s) which shows a feedback bandwidth of 100 kHz. 



Ideally, the Hclmholtz coil current supply would intro- 
duce no additional frequency dependence to P(s), so that 
P(s) would be determined solely by atomic dynamics. In 
the limit, (F z ) <C |F|, which is maintained in closed-loop, 
this would be [using (F z (t)) = exp(— t/T2) sin(ojit)], 

P(s) w 2 7 VM|F|- (8) 

s 

and proportional feedback [constant C(s)] would pro- 
vide good tracking. However, it is not always possi- 
ble to construct a current supply with sufficient band- 
width. Finite available supply power places an upper 
bound on the closed-loop bandwidth for driving an in- 
ductive load, and a more intelligent C(s) design is neces- 
sary. Our Helmholtz coil supplies were constructed using 
high-power (500 W) operational amplifiers and displayed 
73 degrees of phase delay at 100 kHz. With proportional 
control, this would be insufficient phase margin. 

Fig. shows the measured transfer function, P(s), for 
the combined atomic ensemble, Faraday polarimeter and 
feedback coil system. It was generated by utilizing a net- 
work analyzer to perform a swept sine analysis of P(s). 
For each data point, the analyzer was triggered follow- 
ing preparation of a coherent spin state according to the 
cooling and optical pumping procedure described above. 
Although P(s) is approximately an integrator, it displays 
substantially larger phase delay at higher frequency; a fit 
of the transfer function yields the model, 

= 1.6xlQ4(8.0xlQ5- s ) 
K ' s 2 + 4.1 x 10 5 s + 4.0 x 10 9 ' K ' 

By factoring this transfer function into its minimum- 
phase and all-pass components, P(s) = P mp (s)P ap (s), 
a stabilizing feedback controller was obtained using Eq. 
(0). W(s) was chosen to be a single-pole (Butterworth) 



low-pass filter with / c = 1 MHz which yielded the trans- 
fer function, C(s), in Fig. |^A). The feedback controller 
was implemented using high bandwidth analog electron- 
ics and resulted in the closed-loop transfer function, T(s), 
depicted in Fig. [JB). The significance of T(s) is that 
|| 1 — T(s) | provides a measure of the tracking error, and 
thus the error in the parameter estimation. 

Fig.^shows a demonstration of the closed-loop param- 
eter estimation procedure for both stationary and time- 
dependent magnetic fields generated by an auxiliary y- 
axis Helmholtz coil (refer to Fig. |2J). Plot (A) indicates 
the effect of real-time feedback on the polarimetry pho- 
tocurrent, y(t). At time, t — 0.5 ms, a field of B y = 50 
mG was turned on, with feedback disabled, resulting in 
atomic Larmor precession. However, when the feedback 
loop was closed at t = 1 ms, it acted to null out the ap- 
plied field and lock the Bloch vector into the xj/-plane. 
This simultaneously provided an estimation of the ap- 
plied field, b = b c = —0u(t)y. Plot (B) demonstrates 
our capability to track time- varying fields within the feed- 
back bandwidth of 100 kHz; here, feedback was enabled 
the entire time during a 5 kHz bandwidth applied field. 

To demonstrate robustness, the parameter estimation 
error was measured as a function of atom number, N, by 
varying the MOT loading time, t^, according to a cali- 
bration of N versus ti,, (accurate to 40%) obtained from 
resonant fluorescence imaging. The estimation error, Ab, 
was computed from ensemble averages, i?[||b(t) — b(i)|| 2 ], 
over 100 replicate measurements for each sampled value 
of according to the semi-norm, 

\\h(t)\\ 2 = ±£\b(t)\ 2 dt. (10) 

As seen in Fig. 03 Ab was essentially unaffected by fluc- 
tuations in A^ spanning three orders of magnitude. The 
residual estimation uncertainty of approximately 713 /iG 
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FIG. 4: Single-shot, closed-loop parameter estimation of a stationary (A) and time-varying (B) applied magnetic field h(t). 
In (A) feedback is enabled at t = 1 ms (dotted line). The labels identify several artifacts in the data: (1) cross-talk between 
the driving and feedback coils, (2) background field fluctuations due to power supply noise at 51.3 kHz, and (3) noise in the 
driving field that is also revealed by the estimator (4). 



is the result of tracking error due to finite controller gain 
which provides motivation for higher bandwidth feedback 
systems. A closed-loop estimator with a 1 MHz unity- 
gain point is expected to provide ~10 nG field resolution 
similar to current magnetometers |l5) . 

Although the current experiment did not produce 
an appreciable degree of spin-squeezing, it has demon- 
strated a connection between feedback and robustness 
(to inevitable atom-number fluctuations) that will apply 
equall y w ell in future experiments with improved sensi- 
tivity jlfj. Hence the closed-loop methodology we advo- 
cate should enable — and may even be essential for — 
the implementation of proposals to utilize conditional 
spin-squeezing for sensitivity beyond the Standard Quan- 
tum Limit in various applications of spin resonance with 
cold atoms, such as magnetometry atomic frequency 
standards, and matter-wave gravimetry. 
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FIG. 5: Robustness of the single-shot, closed- loop estimation 
to atom number fluctuations (refer to text). 
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